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By characterizing the dynamics of idealized lattice models with a tunable kinetic constraint, we 
explore the different ways in which dynamical facilitation manifests itself within the local dynamics 
of glassy materials. Dynamical facilitation is characterized both by a mobility transfer function, 
the propensity for highly-mobile regions to arise near regions that were previously mobile, and by a 
facilitation volume, the effect of an initial dynamical event on subsequent dynamics within a region 
surrounding it. Sustained bursts of dynamical activity - avalanches - are shown to occur in kineti- 
cally constrained models, but, contrary to recent claims, we find that the decreasing spatiotemporal 
extent of avalanches with increased supercooling previously observed in granular experiments does 
not imply diminishing facilitation. Viewed within the context of existing simulation and experimen- 
tal evidence, our findings show that dynamical facilitation plays a significant role in the dynamics of 
systems investigated over the range of state points accessible to molecular simulations and granular 
experiments. 



I. INTRODUCTION 

When glassy materials, such as supercooled liquids, 
dense colloidal suspensions or driven granular materi- 
als, are cooled or compressed towards the glass or jam- 
ming transitions, the motions of their constituent parti- 
cles become increasingly correlated in space and time p]- 
[3] . This phenomenon, known as dynamical heterogeneity 
(DH) , is a universal property of glassy materials [4-6] and 
is thought to have a direct connection with the anoma- 
lous transport properties of these systems near the glass 
and jamming transitions [7HT4]. 

Despite its ubiquity, the microscopic mechanism of DH 
remains uncertain. Recent insights into this problem ob- 
tained through molecular simulations and granular ex- 
periments have demonstrated that DH can be decom- 
posed into smaller dynamical sub-units. Gebremichael, 
et al showed that strings of highly-mobile particles [15] 
are made up of shorter microstrings whose character does 
not change with temperature, although the strings them- 
selves grow with supercooling p~6j[T7]. Candelier, Dau- 
chot, Biroli and co-workers showed that clusters of par- 
ticles undergoing nearly-simulataneous cage escapes co- 
alesce into larger mobile clusters [18-20 ]. A comprehen- 
sive study of elementary excitations [21] showed that DH 
builds up from localized excitation dynamics involving 
the collective displacements of only a handful of parti- 
cles spanning a few molecular diameters. These observa- 
tions imply the existence of spatiotemporal correlations 
between dynamical subunits, where dynamical events fa- 
cilitate subsequent dynamics nearby in space, giving rise 
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to large-scale DH over time [9j [10] . The degree to which 
such dynamical facilitation (DF) plays a role in the relax- 
ation mechanism of glassy materials is disputed, even for 
systems for which the local dynamics can be resolved di- 
rectly. This lack of agreement stems from uncertainty re- 
garding how DF manifests itself on the microscopic level, 
and as a result, several different methods have been pro- 
posed for measuring DF, leading to contrasting interpre- 
tations. In particular, several studies [2TH23] report that 
DF is present at all supercooled state points, whereas 
others argue that DF must be augmented by another 
structural relaxation mechanism [19] or is only relevant 
over a narrow range of state points [20] . 

Here, we apply several dynamical characterization 
schemes, originally proposed in the context of measur- 
ing DF in molecular simulations and granular experi- 
ments, to kinetically constrained lattice models [24], for 
which DF is the primary relaxation mechanism by con- 
struction. These quantities are also measured for hy- 
brid models that range between a hard dynamical con- 
straint (pure facilitation) and a non-interacting lattice 
gas, which allows us to study the effect of delocalized soft 
rearrangements in violation of DF of the type proposed 
by Ref. [19] . To allow for comparison with molecular sim- 
ulations and granular materials, we formulate a displace- 
ment field for kinetically constrained models that is anal- 
ogous to coarse-grained particle displacements within a 
small region of space for particulate systems [21 . We 
show that mobility transfer correlations, based on the ex- 
change of mobility amongst neighboring spatial regions, 
and facilitation volumes, based on correlations between 
dynamical events and subsequent dynamics within the 
surrounding sub- volume, accurately reflect the degree of 
DF within a given system. We verify that avalanches 
of sustained dynamical activity arise from facilitated dy- 
namics, as predicted by Refs. [18-20 , but in contrast 
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with the interpretations of Ref [20 , we find that the ten- 
dency for avalanches to contain fewer dynamical subunits 
with increased supercooling does not imply diminishing 
facilitation. We predict that a decreasing role of DF 
could be detected by mobility transfer correlations that 
go through a maximum, or facilitation volumes that ex- 
hibit little growth with increased supercooling, although 
these behaviors have not been observed in either simula- 
tions or experiments. 



II. MODELS AND SIMULATIONS 

DF presumes that the glassy material contains local- 
ized soft-spots, or excitations, that allow for local struc- 
tural rearrangement. These rearrangements facilitate the 
birth and death of excitations nearby in space, thereby 
facilitating nearby motion at a later time. This physical 
picture is encoded into a class of kinetically constrained 
models (KCMs) that represent excitations as spins in 
a non-interacting lattice gas, where lattice sites change 
state only in the presence of neighboring excitations. As 
a result of this dynamical constraint, these simple mod- 
els exhibit the complex structural relaxation behavior of 
glassy materials [10] [25j [26] , although their equilibrium 
thermodynamic behavior is trivial. 

We consider one-dimensional (d=l) systems with TV = 
1024 lattice sites occupying one of two states rti = {0, 1}, 
where 1 (0) represents an excited (unexcited) state. The 
system Hamiltonian is given by H = and thus 

excitations are present with equilibrium concentration 
c = (rti) = (1 + e^) _1 . The inverse temperature is given 
by /3 = l/k&T with Boltzmann's constant &b taken as 
unity. We consider systems with a directional dynamical 
constraint defined by the Id east model [24] [27], where 
sites with rii = 1 can facilitate the adjacent site r^+i. 
Structural relaxation in the east model follows from hi- 
erarchical dynamics of the type proposed by Ref. [28] 
and reported in simulations of atomistic supercooled liq- 
uids [21] . Below the glassy dynamics onset temperature, 
T Q , this structural relaxation law is in good agreement 
with simulation and experiment [25]. We approximate 
T Q here as the maximum T for which the systems ex- 
hibits significant four-point correlations, T Q « 1. The 
structural relaxation time r is defined as the mean time 
required to for relaxed regions to span the mean distance 
between excitations, £ = c~ d . 

In addition to facilitated moves, we allow for moves 
that violate the kinetic constraint with probability 
exp(— f)U so ft). The parameter U so it represents the energy 
barrier for soft, delocalized dynamics, and is systemati- 
cally varied throughout the study. Taking into account 
the constraint and detailed balance, the transition rates 
for a given lattice site are given by fco->i = e ~^( n i-i + 
exp(-/3Z7 so f t )) and k^o = + exp(-/3Z7 so f t ). In ad- 
dition to the directional east model, in many cases, we 
have compared our results with the non-directional FA 
model [29], and verified that they are qualitatively simi- 



lar, although these results are not presented here. 

III. DYNAMICS OF KINETICALLY 
CONSTRAINED MODELS 

DF is trivial to measure in KCMs because the exci- 
tations themselves are directly observable. Such direct 
measurement is not currently possible in molecular sim- 
ulations or related experiments, because the precursors 
to excitation dynamics, if they exist, are not yet known. 
Instead, DF must be inferred from dynamical quanti- 
ties. Dynamics in molecular systems correspond to lo- 
cal particle rearrangements, which, in turn, correspond 
to rearrangements in the underlying positions of excita- 
tions. Thus, by analogy, dynamics in KCMs correspond 
to changes in microstate, or kinks. A kink occurs at site 
i and time t if ^(t) = rii(t) — rii(t — St) is ±1, where St is 
an elementary time step. Although excitations connect 
throughout space and time, kinks become disconnected 
when viewed on short time scales, as illustrated by the 
main panel and inset of Fig. for a long east model 
trajectory. 

Particles in molecular systems surge back and forth be- 
fore eventually sticking to new positions and giving rise 
to long-lasting displacements [21] [30] . Similar behavior 
is observed for the east model, and other hierarchical 
KCMs [10], where the majority of kinks are quickly re- 
versed, producing only fleeting changes in the arrange- 
ment of excitations. This behavior can be observed in 
the inset of Fig. [ljt, which shows a short excerpt from 
a longer trajectory. Ref. [21 introduces a method for 
distinguishing between these fleeting reversible motions 
and excitation dynamics in molecular systems by selec- 
tively sampling events involving particle displacements 
that are maintained for a significant period of time [21]. 
Analogously, for KCMs, we consider a subset of enduring 
kinks for which a change in occupancy is maintained for 
a sojourn time, t s . Here, we take the sojourn time as the 
mean time scale for dynamical exchange events, (r x ) [31], 
which, for atomistic models, is on the order of the plateau 
times used to characterize excitation dynamics [21] [32] . 

An enduring kink at site i and time t is indicated by a 
function hi(t), defined below in terms of three binary op- 
erators. These operators are described in the color-coded 
schematic shown in Fig. The first operator, At), 
requires that kinks on both ends of the trajectory endure 
longer than a sojourn time, t Sl 

Si(t,At)= Y[ ^(t + sAt/2) 
-={-i,i} 

(1) 

x Y[ [1 - K 2 i(t + sAt/2 + st')] . 

t'=5t 

The first product requires that kinks occur at both ends 
of the enduring kink event, spanning from t — At/ 2 to 
t + At/ 2. The second product ensures that no kinks 
occur within a time t s prior to the first kink or within a 
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FIG. 1: Dynamics of kinetically constrained models, (a) The panel depicts excitations and kinks for a long trajectory in 
the east model. While the excitations connect throughout space and time, the kinks are temporally separated. The inset 
illustrates that this feature becomes apparent on shorter time scales, (b) The panel depicts a short trajectory for the east 
model. Enduring kinks (long-lived changes in the arrangement of excitations) are depicted by red x's. At the bottom of the 
panel, lattice sites are colored according to /i;(0,£), the number of enduring kinks at site i up to the time t. (c) The schematic 
depicts a small region of space-time, containing cells that span an elementary time step St, with dark shaded cells containing 
excitations, rij{t) — 1, and white cells containing no excitations, rij(t) — 0. Black cells highlight excitations at a lattice site z, 
which exhibits an enduring kink at time t, denoted by a red x. Enduring kinks are defined by the binary indicator function 
hi(t), which itself is defined in terms of three binary operators. These operators, <Si(£, At), T(t, At), and (D(t, At), are described 
pictorially by arrows and text in the panel, and color-coded blue, red, and green, respectively. Corresponding mathematical 
definitions are provided in the text. 



time t s after the final kink. When these conditions are 
satisfied, <Si(t, At) = 1. Otherwise, it equals zero. 

The transient portion of the trajectory is defined such 
that at least one kink occurs within every sliding time 
window of size t s — St between time t — At/2 and t + At/2, 



7^, At) 



t-|-At/2-t s 

n 

t'=t-At/2 



t'+t s 

1-6 \ 

t"=t'+5t 



(2) 



When this criterion is satisfied, 7i(t, At) = 1. Otherwise, 
it equals zero. The transient portion of the trajectory 
must contain an odd number of kinks, such that the event 



results in an overall change in state, 

(*+At/2-t B 
t'=t-At/2 



(3) 



Oi(t, At) = 1 when an odd number of kinks have occured 
between time t — At/2 and t + At/2 and zero otherwise. 
The path functional for an enduring kink is then given 

by 

hi(t) = s i(t, At)75(t, At)Oi(t, At). (4) 

At 

The term in the summation is a product over all of the bi- 
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nary operators defined above. The summation is carried 
out over all possible durations, At, where, by construc- 
tion, only one value of At can satisfy all of the operators 
simultaneously. The function hi(t) is therefore itself a 
binary operator and equals unity if and only if all three 
conditions are satisfied. Otherwise, hi{t) equals zero. 

The displacement at a lattice site i over a time window 
t r — t is approximated as a simple sum of enduring kinks, 



(5) 



In molecular systems, the total displacement of a given 
particle does not typically scale linearly with the num- 
ber of discrete displacements, since it is unusual for dis- 
placements to occur in exactly the same direction. Thus, 
the mobility field may be more accurately described by 
[ii(t,t f ) p , where the power p depends on the fractal di- 
mensionality of particle diffusion [33 , 34 j. For simplicity, 
we assume that p = 1, as this does not effect the qualita- 
tive behavior of the measurements performed here. More 
realistic mappings might also account for temporary dis- 
placements that arise from reversible surging. We ob- 
serve that quantities based on jHi(t^t f ) exhibit artifacts 
for t < t s , due to coarse-graining on a time scale t s . This 
does not affect the measures considered here, which in- 
volve significantly longer time scales. 
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IV. MEASURING DYNAMIC FACILITATION 

With the dynamics of KCMs defined, we study the 
properties of three different measures of DF, originally 
proposed in context of simulations of glass-forming liq- 
uids and granular materials in experiment. The measure- 
ments are presented chronologically, and compared based 
on their proficiency at detecting DF, their temperature 
variation, and their ability to distinguish between sys- 
tems with a hard dynamical constraint and those with 
softened dynamics. 

A. Mobility Transfer Function 

Vogel, Glotzer and coworkers have proposed a mobil- 
ity transfer function based on the probability of observing 
new mobile particles near particles that were previously 
mobile. For facilitated dynamics, this probability exceeds 
the uncorrelated case, and this trend becomes more ap- 
parent with supercooling, as c decreases [22j [23] . 

To measure this function for KCMs, we begin by intro- 
ducing a small perturbation to the mobility field at each 
lattice site, 

Mi(t' ,t") = Hi(t' ,t") + W(0,e). (6) 

^/(0, e) is a uniform random number on the interval [0, e], 
where e is a small number <C 1. This perturbation en- 
sures that no two lattice sites exhibit exactly the same 



FIG. 2: Mobility transfer function, (a) The main panel 
shows a representative curve for Pf(R), the probability of 
observing a newly-mobile site a minimum distance R away 
from a previously- mobile site. The red dashed curve shows 
the uncorrelated result, obtained by choosing the initial "mo- 
bile" sites at random. The inset shows the value of F(t), the 
integrated difference between the curves for each t, which rep- 
resents the degree of DF. (b) Comparison of the peak value 
of F(t), ^? max (7 1 ) 5 as a function of temperature for different 
values of the softness parameter, U so ft 

mobility, which becomes important in certain limits of 
the analysis. We define mobile sites according to the fol- 
lowing binary operator 

mi(t',t") = 6 [Mi(t',t") - Mcut(i',0] • (7) 

Here, 6, is the Heaviside step function, and /i cut is chosen 
to include the 5% most mobile sites during a time window 
spanning from t' to t" 

/•Meut(t , ,0 N 

/ Yl 6 [ Mi ~ ^ d »' = i°' 9 $ N \ ' ( 8 ) 

Here, |_- • • J denotes the floor function. Sites with 
rrii(t' \t") = 1 comprise the highly-mobile subset. For all 
other sites, rrii(t f : t ff ) = 0. The random perturbation en- 
sures that the 5% most-mobile sites are well-defined, even 
when several sites exhibit the same value of [ii(t',t"). A 
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set of newly-mobile sites is denned by 

Wi(1?X) = [l-m i (2t / -t",t')}m i (t',t"), (9) 

where the interval is symmetric about t' . Here, 
Wi(t' ,t") = 1 for sites that are newly mobile, and zero 
otherwise. For each newly-mobile site, we measure the 
minimal distance P to a previously- mobile site. That is, 
we find Ri = \D(i,j)\ that minimizes 

Ri = arg min ( - 1 ... rnj {t' , t"K(*' , t") \ , (10) 

\D(i,j)\ I \ V \hJ)\ ) 

where, D(i,j) is the distance between sites, 

D = \i-j\-NQ(\i-j\- N/2) . (11) 

Here, the second term is needed due to periodic bound- 
ary conditions. The probability distribution of R as a 
function of the time window t = t" — t 1 ', Pf(P, t), is com- 
pared to the uncorrelated distribution, Pp(P, t). For this 
distribution, Mi{t' \t") is replaced by U(0, e), such that 
the N — [0.95N\ "mobile" sites are chosen at random. 

Both P F (P, t) and Pp(P, t) are plotted in Fig. [2^ for a 
typical state point and value of t. The value of Pf(P, t) 
exceeds Pp(P, t) for small values of P, indicating a pref- 
erence for newly-mobile sites to arise close to previously- 
mobile sites. This preference is further quantified by the 
mobility transfer function, 

F{t)=<TP ¥ {R',t) I f P^,t). (12) 

We use P cut = 2 lattice sites, but we obtain similar re- 
sults for < P cut < 4. The inset of Fig. [2^i shows the 
time evolution of F(t) as a function of temperature. We 
find that, as for atomistic and molecular glass formers, 
F(t) exhibits a peak value, F ma *(T), that shifts to later 
times and grows with decreasing temperature. Fig. [2^i 
shows that F max (T) is roughly proportional to the mean 
distance between excitations £ for the east model. This 
indicates that values of P become increasingly dominated 
by facilitation as c decreases, and the probability of ran- 
domly encountering nearby mobile sites in the absence of 
facilitation decreases proportionally. 

We find that this relationship breaks down when soft- 
ness is added to the system. In particular, Fig. [2}d shows 
that softened systems exhibit a crossover, where F max (T) 
first increases and then decreases as a function of tem- 
perature. This crossover in F max (T) corresponds to a 
crossover in the relaxation mechanism, where soft delo- 
calized dynamics - given an Arrhenius temperature de- 
pendence for our models - becomes more probable than 
DF, which has a super-Arrhenius temperature depen- 
dence. Thus, the temperature variation of F ma *(T) gives 
a qualitative measure of softness, provided that soft de- 
localized relaxation becomes dominant at lower T. The 
fact that atomistic and molecular simulated glass formers 
do not exhibit such a crossover [22] [23] indicates that soft 



delocalized relaxation does not dominate the dynamics, 
at least over the range of temperatures studied. We find 
that finite size effects, which appear when £ approaches 
the order of the system size, give rise to a qualitatively 
similar crossover behavior, and therefore care should be 
taken when interpreting simulation or experimental re- 
sults in future studies. 



B. Avalanches 

Candelier, Dauchot, Biroli and coworkers [T8H20] 
have extensively characterized the dynamics of so-called 
"cage-escapes," defined by particles that obtain a new 
center of vibrational motion. Cage escapes coalesce into 
clusters in space and time that resemble the excitation 
dynamics described in Ref. [21]. The tendency for clus- 
ters of cage escapes to facilitate one another is character- 
ized by a distribution of waiting times between adjacent 
clusters, t\. For KCMs, n is measured according to the 
following argument minimization over time t: 

n = argmin{-(l/t)(l-^(0)/€ nbr (t)])}. (13) 

t 

Here, nbr stands for neighboring sites, j with \D(i,j)\ = 
1. In simulated supercooled liquids and granular ma- 
terials in experiment, the probability distribution of ti, 
Pi(t), resembles the superposition of two exponential dis- 
tributions [18]. The time constant r corr characterizes 
short-time exponential behavior. It is speculated that 
this duality of t\ distributions implies two distinct phys- 
ical mechanisms, with short lag times arising from fa- 
cilitated dynamics and long lag times arising from soft 
delocalized dynamics [19] [20]. This argument implies 
that the east model, or any other purely-facilitated KCM, 
should exhibit single-exponential behavior in P\(t\ since 
DF is the only relaxation mechanism in these systems. 
However, as illustrated in Fig. (3^i, the distribution of 
Ti for the east model resembles the double-exponential 
behavior observed for model liquids and granular mate- 
rials. Moreover, the inset of Fig. [3^i shows that intro- 
ducing soft, delocalized dynamics tends to diminish the 
two-exponential effect, rather than amplify it. 

Our findings for the east model highlight the impor- 
tant distinction between excitations and their dynamics. 
While excitations necessarily connect throughout space 
and time for facilitated models, their dynamics becomes 
intermittent at lower T, giving rise to a wide range of 
lag times between events. This can be rationalized in 
terms of persistence and exchange times, quantities that 
arise in the context of transport decoupling in super- 
cooled liquids [31] [32]. In KCMs, persistence times, r p , 
are the time scales over which randomly-chosen lattice 
sites exhibit their first kink. Exchange times, r x are the 
lag times between subsequent kinks at the same lattice 
site [31]. Like the exchange and persistence time distri- 
butions, P x (t) and P p (t), the distribution of facilitation 
lag-times Pi(t) exhibits rich behavior, and is described 
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FIG. 3: Facilitation waiting time distributions (a) Proba- 
bility distribution of lag times between adjacent kinks, Pi(t), 
as a function of temperature for the east model. The inset 
shows the effect of including soft, non-localized relaxation for 
one temperature. The parameter P Q corresponds to the value 
of the first bin of the Pi(t) histogram, (b) Exchange and per- 
sistence time distributions P x (t) and P p (t) for the east model 
at two different temperatures compared to Pi(t). 



only approximately by the sum of two exponentials, at 
least for KCMs. This is demonstrated by plotting P\(t) 
on a log-linear scale, as depicted in Fig. ^jp. The fig- 
ure shows that P\{t) involves a combination of persis- 
tence and exchange-like time scales. As P*(t) and P p (t) 
decouple at low temperatures [31j [32], P\{t) becomes 
more separated between long and short lag times. This 
may explain the tendency for clusters of cage escapes to 
form large avalanches at higher temperatures, but be- 
come more sporadic at low temperatures, when clusters 
are grouped into avalanches based on the short-time ex- 
ponential time scale of Pi(t), as described in Ref. [20] . 

We explore this possibility further by explicitly defin- 
ing avalanches for the east model. In the analysis de- 
scribed above, we considered waiting times for kinks K>i(t) 
to place emphasis on the fact that the observed double- 
exponential behavior in P\(t) is robust, even for unpro- 
cessed dynamics. In the analysis that follows, we con- 
sider what we believe to be a more realistic mapping 
of KCMs to clusters of cage escapes by replacing kinks 
Ki(t) with enduring kinks hi(t) in Eq. (|13|. This does 



not change the qualitative behavior of Pi(t), but rather 
truncates the distribution for very short lags, less than 
the sojourn time, t\ < t s . Avalanches are defined ac- 
cording to Refs [18-20 by grouping enduring kinks that 
are adjacent in both space and time. Space- time is dis- 
cretized into points (i, t) with temporal lattice spacing St. 
Points (i,t) and (j,t f ) belong to the same avalanche a if 
they contain enduring kinks separated by \t — t'\ < r corr 
and \D{i,j)\ < r corr . Each point belongs to exactly one 
avalanche by construction. The set of all avalanches a for 
a state point is denoted A. The cutoff r corr is taken to be 
two lattice sites and r cor r is set by the short-time expo- 
nential cutoff of P\(t) for enduring kinks. Because P\{t) 
is only approximately exponential, as described above, 
Tcorr depends somewhat on the histogram bin size, Ar±. 
For Ati sufficiently small (Ati ~ 10) and chosen con- 
sistently for all data sets, the qualitative behavior of 
avalanches is not affected. Fig. [4^i compares r cor r to the 
average time duration of an avalanche, Tf ac , defined by 



max 

(i,t)ea , (j,t f )ea uv 



i-(j,t f )'i\} 



(14) 



Here, ('")a denotes an ensemble average over 
avalanches, a, and t is a unit vector (0, 1) that projects 
each point into the temporal dimension. For compari- 
son, Fig. also shows Ti/ 2 , the average minimum time 



required 



or 1/2 the lattice sites to exhibit a kink, 



arg mm 




<N/2 . (15) 



This time scale is similar to r a for particulate sys- 
tems [20]. While r CO rr does not vary much with temper- 
ature, Tf ac decreases as T is lowered, in agreement with 
the results obtained for the air-driven granular system 
studied in Ref. [20]. 

The sequence of trajectories in Fig. illustrates the 
qualitative evolution of avalanches with changing tem- 
perature. At high T, dynamics occurs within a single 
avalanche. As temperature is lowered, avalanches become 
increasingly intermittent and spatially separated. This 
is quantified by measuring the probability of observing a 
spanning avalanche with Tf ac > Ti/ 2 as a function of T. 
As shown in the inset of Fig. [4^i, the spanning probability 
crosses over sharply at around T = 0.6. This crossover 
behavior is robust, although we find that the position of 
the crossover depends on the parameters chosen to de- 
fine the avalanches. This indicates that a decreasing spa- 
tiotemporal extent of avalanches with supercooling does 
not preclude DF but rather seems to represent an intrin- 
sic behavior of purely-facilitated systems. 

While the temperature variation of avalanche time 
scales closely mirrors that of molecular simulations and 
granular materials, the variation of spatial length scales 
is more difficult to compare. Close inspection of the low- 
est temperature state point depicted in Fig (4}d reveals 
that the average spatial extent of the avalanches becomes 
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FIG. 4: Avalanches in KCMs (a) Time scales r CO rr, Tf ac and Ti/ 2 as a function of temperature for enduring kinks in the 
east model (see text for definitions). The inset shows the probability of observing a spanning avalanche of time duration T1/2 
or greater as a function of temperature, (b) Avalanches in east model at three different temperatures. Enduring kinks are 
depicted by circles and colored according to the avalanche to which they belong. Trajectories span approximately 3ti/ 2 . (c) 
The four point correlation length £4 and the avalanche spatial length scale £ ava as a function of temperature. 



smaller than the typical dynamical length scale £ (the dis- 
tance between excitations, shown in grey). This is quan- 
tified by computing the mean end-to-end distance of an 
avalanche [20] , 

max {\(i,t) ■ s - (j,t') ■ s\}\ . (16) 

(t,t)£a , (j,t')ea I A 

Here, s = (1,0) is a unit vector that projects each point 
into space. The avalanche length scale is compared to the 
dynamical correlation length £4 [35, 36 , defined in terms 
of 



X4(t)=N (q(t)) 2 - (q(tf) 



(17) 



is zero 



The indicator function qi(t) — S 

if at least one kink has occurred at lattice site i over a 
time window t and one otherwise. We obtain £4 from the 
peak value of X4(t), which depends on temperature, £4 = 
j-^max^)] 1 /^ yy e opt to define the four-point length scale 
dynamically for closer correspondence with particulate 
systems; however, t is directly proportional to £4 and 
thus the quantities are interchangeable. 

At high temperatures, for which avalanches span the 
system, we find that £ ava exceeds £4. This behavior is 
also observed for the granular system studied in Ref. [20] . 
There, it is argued that large avalanches on length scales 
exceeding £4 result from a collection of dynamically inde- 
pendent events occurring on a scale £ ava ~ £4. Following 
this argument, we estimate £ ava « £4 in this range. Be- 
low a crossover temperature, £ ava becomes smaller than 
£4 and decreases with supercooling. The low dimension- 
ality of the model studied here may accentuate this trend, 
but this qualitative behavior should hold in general. This 
stands in contrast with the interpretations of Ref [20 , 
where it is speculated that avalanches span a typical dy- 
namical correlation length at all state points, and thus 



represent dynamically-independent events. The available 
data for particulate systems seems to neither prove nor 
disprove this hypothesis. Ref [20 finds that £ ava ~ £4 
near the crossover, but this is true by construction, given 
that a crossover exists and the temperature variation of 
£ ava is not strong, as is the the case for the east model. 
Evidence of dynamically-independent avalanches would 
involve £ ava that grows in proportion to £4 over a range 
of supercooled state points, but, for the narrow range of 
state points available for particulate systems, £ ava neither 
definitely grows nor definitively shrinks. We thus leave 
our findings as predictions to be investigated through fu- 
ture study. 



By studying avalanches in KCMs, we arrive at a some- 
what different interpretation regarding the specific re- 
lationship between avalanches and DF than that previ- 
ously articulated. In particular, we find that the ob- 
served phenomenology previously thought to be incon- 
sistent with DF - long waiting times between avalanches 
and the decreasing spatiotemporal size of avalanches with 
supercooling - follows from facilitated models and does 
not imply either long-ranged correlations between exci- 
tations [19] or a breakdown of DF at low temperatures 
or high packing fractions [20] . Despite these differences, 
our interpretations are in agreement with the overarching 
spirit of Refs. [T8H20] - that is, we find that avalanches 
do represent a type of facilitated dynamics and that lo- 
calized dynamics in particulate systems do not map to 
directly to excitations in KCMs. Our findings reinforce 
the idea that avalanches are a robust phenomenon war- 
ranting future study, particularly with regards to their 
dependence on temperature and dimensionality, as well 
as their dynamic independence and connection to excita- 
tion dynamics. 
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C. Facilitation Volume 

The mobility transfer function and avalanches de- 
scribed in the previous sections involve the practical dif- 
ficulty of defining cutoff values that, if chosen incorrectly, 
can affect the qualitative outcome of the measurements. 
Ref. [21] introduces an alternate, parameter- free mea- 
sure of DF known as the facilitation volume, which quan- 
tifies the overall impact of an initial rearrangement on 
the subsequent mobility field [21 . The time-dependent 
displacement field conditioned on an enduring kink at a 
tagged site 1 is given by, 



T = 0.48 



N 



£(r,t,t') 



MO)Eft(M'Wr-Wj') 
<M0)>n(r) 



-. (18) 



Here r is an integer in the range [0, N/2] and n(r) is the 
probability density of observing a distance r between lat- 
tice sites, n(r) = [l + 0(r)]. We consider the behavior 
of /i(r, t, t') over a time window spanning from At/ 2 to 
t, where At/2 demarcates the completion of the endur- 
ing kink at the origin. Thus, the initial dynamics does 
not factor into the value of the displacement field. The 
quantity /i(r, t, t') is closely related to the distinct contri- 
bution to the four-point susceptibilities %4 and xu, often 
used to characterize DH [35j|37]; however, /i(r, t, t f ) fo- 
cuses specifically on dynamical correlations with initial 
dynamics, whereas four-point functions are sensitive to 
large-scale dynamical correlations that build up from the 
initial dynamics over time. That is, /i(r, t,t') is a three- 
point correlation function [38 , rather than a four-point 
correlation function. 

The behavior of /i(r, At/2, t) is plotted in Fig. |5^i as a 
function of the time window t for the east model. The in- 
set shows that, for t s < t < 5r, /i(r, At/2, t) peaks near 
position of the previously-enduring kink at the origin, 
and decays to (/i(0,t)) as r becomes large. The main 
panel shows the same curves normalized by (/i(0,t)), 
highlighting the excess displacement relative to the un- 
corrected result obtained in absence of the initial dy- 
namics. The curves plateau at r « £/2, indicating that 
mobility correlations span a maximum range £ ~ 1/c. 
The fact that the length scale is relatively constant for a 
large range of time scales reflects the non-zero probabil- 
ity of observing a chain of excitations that spans £, even 
on very short time scales. 

The facilitation volume is defined as the sum of the 
excess mobility over all r, 



v F (t) 



N/2 

£ 

r=0 



fi(r,At/2,t) 
|>(At/2,t)> 



(19) 



The behavior of vp(t) is plotted in Fig. |4Jd as a func- 
tion of temperature. For each temperature, the function 
peaks near the sojourn time t s and decays toward zero 
as t becomes much greater than r. This is similar to the 




10 100 1000 10000 100000 

t 

FIG. 5: Mobility field and facilitation volume, (a) Time 
evolution of the mobility field jl(r, At/2, t) for the east model 
for a single temperature, T — 0.48. For both the main panel 
and inset, the time window t spans from t ~ t s (blue) to 
t « 5r (red), (b) Facilitation volume VF(t) as a function of 
temperature. Temperatures span the same range as in Fig.[2^L 
The inset shows the peak value of ^(t) at each temperature 
i>F ax (T) for different values of softness, controlled by U so ft- 



behavior observed for atomistic models, exception that, 
for atomistic systems, ^f(^) peaks near the a relaxation 
time [21]. This discrepancy may arise due to the non- 
discrete nature of particulate systems, where vibrational 
motions dominate the displacement field at very short 
times. Regardless of the characteristic time chosen in 
the range t < r, vp(t) grows with decreasing tempera- 
ture, and, in particular, scales as £. This is illustrated 
by the inset of Fig. [5]3, which shows the peak of vp(t) 
as a function of temperature, ^p ax (T). Due to statistical 
uncertainty, it is unclear whether the facilitation volumes 
reported for atomistic systems in Ref. [21 scale approxi- 
mately with £ d , but this relationship should be explored 
in more detail in the future. 

Fig. also shows the variation of ^p ax (T) with the 
characteristic energy for soft delocalized relaxation, C/ so ft • 
Allowing for non-facilitated dynamics moves adds back- 
ground noise to the displacement field, which diminishes 
the excess mobility and reduces ax (T) accordingly. For 
large amounts of softness, relaxation occurs mostly in 
the absence of facilitation, and Vp ax (T) no longer ex- 
hibits significant growth with decreasing T. In contrast 
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to F max (T), the value of vj? ax (T) does not exhibit a tem- 
perature crossover for any value of f7 so ft- This is related 
to the fact that v]? ax (T) is normalized by (jx(At/2, t)), 
which decreases as a function of T, whereas i^ max (7") i s 
normalized by the same function for all T. Thus, our 
findings imply that large values of ^p ax (T) that grow 
with supercooling, such as those observed in Ref [21], in- 
dicate the presence of significant DF, but do not rule out 
the possibility of some delocalized relaxation as well. 

V. DISCUSSION 

Our results indicate that the current body of litera- 
ture regarding DF implies that DF is present at all state 
points investigated, and becomes increasingly apparent 
with increased supercooling. It remains possible that DF 
could be superseded by another mechanism at currently- 
inaccessible conditions, but evidence for such a mecha- 
nism has not yet been reported. The fact that structural 
relaxation data for a wide range of experimental condi- 
tions corresponds to functional forms predicted for facili- 
tated dynamics [25] seems to indicate that the relaxation 
mechanism does not change, although other interpreta- 
tions are possible [8]. 

Our results regarding avalanches imply that clusters of 
cage escapes are closely related to excitation dynamics 
and might be applied to study DF in the future, partic- 
ularly for experimental systems, or other systems where 
inherent states or time-coarse graining on a very fine scale 



is not possible. In Ref. [19] soft modes [39] were invoked 
to explain the sudden triggering of avalanches at points 
in space with little prior dynamical activity over long pe- 
riods of time. Although our results imply that such long 
waiting times are a natural consequence of facilitated dy- 
namics, the observed connection with soft modes implies 
that these quantities might provide a generic method for 
detecting excitations in absence of their dynamics [40] . 

In addition to the quantities involving DF explored 
here, the method that we derive for a displacement field 
for KCMs might allow for additional detailed mapping 
between the dynamics of KCMs and molecular systems. 
The quantities reported in Ref [21] seem like a particu- 
larly fruitful avenue for future study. 
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